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width of beam and piezo-actuator, m 
electric charge constant of piezo-actuator, m/v 
distance to neutral axis of composite beam measured from 
its lower edge, m 

Young’s modulus of elasticity of piezo-actuator and beam 
respectively, N/m 2 

Young's modulus of elasticity of element i of the beam, 
N/'m 2 

vector of external forces and moments, N or Nm 
area moment of inertia of actuator, bonding layer and 
beam about the neutral axis of the composite beam 
respectively, m 1 * 

area moment of inertia of the element i, m 1 * 
stiffness matrix of the element i 
overall stiffness matrix of beam-actuator system 
length of element i, m 

external moment acting on i tlri node of beam, Nm 
piezo-electric moment generated by piezo-film, Nm 
number of elements of the beam 

thickness of piezo-actuator, bonding layer and beam 
respectively, m 
deflection criterion 

voltage applied across the piezo-electric film, volts 
external force acting on the i th element of the beam, N 
the linear translation of node i, m 



Greek Letters 


deflection of node i, in or rad. 

8 deflection vector of all nodes of the beam, m or rad. 

ff> piezo-electric strain in piezo-actuator, m/m 

a 1,2,3 bending stresses in piezo-actuator, bonding layer and 
beam respectively, N/m^ 
piezo-electric stress in actuator, N/m^ 



ABSTRACT 


This study deals with the utilization of piezo-electric 
actuators in controlling the static deformation and shape of 
flexible beams. 

An optimum design procedure is presented to enable the 
selection of the optimal location, thickness and excitation 
voltage of the piezo-electric actuators in a way that would 
minimize the deflection of the beam to which these actuators are 
bonded . 

Numerical examples are presented to illustrate the 
application of the developed optimization procedure in minimizing 
structural deformation of beams using ceramic and polymeric 
piezo-electric actuators bonded to the beams with a typical 
bonding agent. 

The obtained results emphasize the importance of the devised 
rational produce in designing beam-actuator systems with minimal 


elastic distortions 


INTRODUCTION 


The construction and operation of large structures which are 
extremely flexible have posed new and challenging problems 
particularly because these structures are intended to provide 
stable bases for observations and communications such as in space 
structures. With strict constraints imposed on the structural 
deformations, it became essential to suppress such deformations 
to minimum through the use of active control systems of one form 
or another. 

Distinct among the presently available active control 
systems are those that rely in their operation on piezo-electric 
actuators. Such systems have proven to be experimentally 
effective in controlling the vibrations of simple structural 
elements such as rectangular beams [1-2] and hollow cylindrical 
masts [33. The effectiveness of these systems is coupled also 
with the light weight, high force and low power consumption 
capabilities of the piezo-electric actuators [4-8]. These 
features rendered this class of actuators to be an attractive 
candidate for controlling structural deformation. 

The present state-of-the-art of this type of actuators has 
been limited to the analysis and testing of their characteristics 
[9-11] as influenced by their geometrical or operational 
conditions. But, no attempt has been made towards selecting their 
optimal geometrical parameters or location which are suited for a 



particular structure subjected to known loading conditions. Such 
a synthesis procedure is essential to the successful integration 
of the actuators into the structure in order to minimize the 
structural deformation. 

It is, therefore, the purpose of this study to devise such 
an optimal design procedure that would enable the selection, on 
rational basis, of the optimum design parameters and location of 
piezo-electric actuators to satisfy certain structural 
deformation requirements. The procedure will take into account 
the effect that the actuators have on changing the elastic and 
inertial properties of the structure to which they are bonded to. 
Furthermore, it considers the effect of physical and geometrical 
properties of the bonding layer on the static characteristic of 
the beam-actuator system. 


THE PIEZO-ELECTRIC ACTUATOR-BEAM SYSTEM 
A. General Layout 

Figure (1) shows a general layout of a flexible beam (A) 
whose deflection is to be controlled by piezo-electric actuator 
(B) bonded to the beam by bonding layer (C). The beam, under 
consideration, can generally be made of several steps which are 
not necessarily of the same thickness or the same material. The 
interfacial nodes between the different steps can be subjected to 



Figure (1) - General layout of the piezo-electric actuator-beam system. 


external forces, moments or both. Further, the degrees of freedom 
of any node can be limited to linear translations, angular 
rotations or restrained completely depending on the nature of 
support at the node under consideration. 

In this study, the beam is assumed to have rectangular cross 
section of constant width b and that its transverse deflection is 
due to the flexural action of the external forces and moments. 

In Figure (1), the piezo-electric actuator B is shown bonded 
to the element i of the flexible beam to form a composite beam. 
When an electric field is applied across the film, then it will 
expand if the field is, for example, along the polarization axis 
of the film and will contract if the two were out of phase. The 
expansion or contraction of the film relative to the beam, by 
virtue of the piezo-electric effect, creates longitudinal bending 
stresses in the composite beam which tend to bend the beam in a 
manner very similar to a bimetallic thermostat. 

With proper selection, placement and control of the 
actuator, it would be possible to generate enough piezo-electric 
bending stresses to counter balance the effect of the exciting 
forces and moments acting on the beam in a way that minimizes its 
structural vibrations. 


B. Model Of An Actuator-Beam Element 


Figure (2) shows a schematic drawing of a piezo-film A 
bonded to an element B of the flexible beam. 

If a voltage v is applied across the film, a piezo-electric 
strain (f is introduced in the film and can be computed from : 

= ( d/ 1 1 )*v ( 1 ) 

where d is the electric charge constant of the film, m/v 

t-j is the thickness of the piezo-electric actuator, m 

This strain results in a longitudinal stress Of given by : 

<7 f = E 1 (d/t 1 )»v (2) 

where is the Young's modulus of elasticity of the film, N/m^ 

This, in turn generates a bending moment M f , around the neutral 
axis of the composite beam, given by : 

— (t<+tp+to— D) 

M f = / o f (b*y)dy (3) 

-(t 2 +t 3 -D) 

where t 2 is the thickness of the bonding layer, m 

t 3 is the thickness of the beam, m 

b is the width of the beam, the bonding layer and the 
piezo-film, m 

In equation (3), D is the distance of the neutral axis from the 
lower edge of the beam which can be determined by considering the 
force balance in the longitudinal direction X of the beam, or : 





dA + JoodA + 
"film bond 
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■'beam 


( 4 ) 


or 




-(t 2 +t 3 -D) 


-(t 3 -D) 


ydy + Epbf ydy + Eobf ydy = 0 


(5) 


-(t 1 +t 2 +t 3 -D) -(t 2 +t 3 -D) 


-(t 3 -D) 


where E 2 is Young's modulus of elasticity of the bonding layer 
E 3 is Young's modulus of elasticity of the beam 

Equation (5) yields the following expression for D : 


D = 


E<| 1 1 ^+E 2 t 2 ^+E 3 t 3 ^+2E-j t^ ( t 2 +t 3 ) +2E 2 t 2 t 3 
2(E<|t‘|+E 2 t 2 +E 3 t 3 ) 


( 6 ) 


Equation (2), (3) and (6) can be combined to determine the 
bending moment Mf generated by the piezo-film on the composite 
beam as follows : 


Mf* = 


d*b*E-|*v*(E2t-|t2+E;3t'|t3+E2t2 +2E2t2t2+E^t ^ ^ 
2* ( E-j t ^ +E 2 t 2 +E 3 t 3 ) 


(7) 


For this composite beam, it can be easily shown [12] that it 
has a flexural rigidity (E^I^) given by : 


E^I^ = ^*1^1 + E 2^2 + e 3 i 3 


( 8 ) 


where I-j, I 2 and I 3 are the area moments of inertia of the film, 
the bonding layer and the beam about the neutral axis 
respectively . 


Let us now assume that the composite beam, shown in Figure 



(2), extends a length 1 between two nodes (i) and (i + IK 
Further, it is assumed that the external forces and V i + 1 as 
well as the external moments M e ^ and M e ^ + i are acting on the beam 
at nodes i and i + 1 respectively. Then, the resulting linear and 
angular deformations of the beam yj. and as well as y i+ -j and 
0 i+1 at the nodes i and i+1, respectively, can be related to the 
loads acting on the element as follows [131 : 
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Equation (9) can be rewritten as : 


F t = K^i CIO) 

where Fj_ is the resultant forces and moments vector acting on 
the beam element i, N 

is the stiffness matrix of the composite beam element 
i, N/m 

®i is the deflection vector of the nodes bounding the 
beam element, m 


Equation (9) constitutes the basic finite element model that 
relates the external loads (V and M g ) and piezo-electric moments 
(M f ) to the deflections (y and 0) of the element as a function of 
its elastic and inertial parameters. 


The equation can be equally used for any element of the beam 
whether it has a piezo-film bonded to it or not. In the latter 
case, Mf is set to zero and flextural rigidity is set equal 
to that of the flexible beam element under consideration. 

C. Model Of The Overall Actuator-Beam System 


q The force-displacement characteristics of the individual 

elements of the beam-actuator system, as given for element i by 
equation (9)> are combined to determine the behavior of the 
^ overall structure. 

The equilibrium conditions of the overall structure will be 
• expressed as : 


External forces and moment forces and moments 

acting on the nodes = > acting on the elements 

• of the overall system ^-^at these nodes 


or 


F 


N+1 N + 1 

I F i : S K t 5i = K 5 

i=1 i=1 


( 11 ) 


where K is the overall stiffness matrix of the system (2n*2n) 


Bathe and Wilson [14], Yang [15] and Fenner [16], for 
example, show how to generate the overall matrix K from the 
stiffness matrices Kj. of the individual elements. 


In equation (11), the external force vector F and the 



displacement vector S are given by: 


F t = [V, M-j V 2 M 2 .. .. V N+1 M n+1 ] (12) 

and 

§ T = [Y 1 Y 2 0 2 .. .. Y N+1 #n+1 ^ ^3) 

where superscript T denotes the transpose of the vector 

D. Solution Algorithm 

For a particular flexible beam configuration, external 
loading conditions and end conditions, the effect of applying 
certain voltage v on a piezo-electric actuator of a specific 
thickness t^ placed at a particular location i on modifying the 
elastic deflection 8 of the beam can be determined by solving 
equation (11). 

In equation (11), the overall stiffness matrix K is modified 
to account for the end conditions and type of restraints imposed 
on the different nodes. If, for example, a node (j) is simply- 
supported and is therefore restrained in the Y direction, i.e. 
yjsO, then the off diagonal elements of 2*j row of the overall 
stiffness matrix K are set equal to zero together with the 
corresponding externally applied force. If the node is also 

restrained in its angular motion then 0j=O and the off diagonal 
elements of the (2*j+1) row of the K matrix are eliminated along 
with the corresponding externally applied moment. 



After such a modification process., the elastic deflection 
of the beam can be computed from : 

5 = K""* *F (14) 

where K**^ is the inverse of the overall stiffness matrix 

Figure (3) shows a flowchart of the solution algorithm. 


OPTIMIZATION OF THE ACTUATOR-BEAM SYSTEM 

The presented analysis algorithm of the actuator-beam system 
is utilized as a basic for the development of the optimization 
procedure for selecting the thickness t-j, excitation voltage v 
and location i of the actuator in order to minimize the 
deflection of the beam. 


The optimum design problem is formulated as follows : 

Find thickness t-j, voltage v and location i 

N+1 

to Minimize U = ( I ( + ^ i 2 ))*bt 1 

i = 1 

Vi = vV 


Such that 



Input 

Mo. of elements of beam N 
Young's Modulus of element E^ 
Moment of Inertial of element 1^ 
Element length 1 


Input 

External forces on elements 
External moments on elements M ei 


Input 

Restraints on each node 


Input 

Properties of actuator 
Young's Modulus E^ 
Charge constant d 


Input 

Actuator's thickness t^ 
Applied voltage v 
Location i 


Compute 

Individual stiffness matrix 
and 

Overall stiffness matrix 


Modify 

External moment vector by 
piezo-electric moments M^ C Eq .(7)3 


Modify 

The stiffness matrix to 
account for restraints 


Compute 

The deflection of nodes C Eq .(14)] 



Figure (3) - Flowchart of the solution algorithm 







v min ^ v ^ v max * 

Y . < Y 
A x N i max’ 

and h < ^max 

The considered objective function aims at minimizing the 
linear and angular deflections of the beam as well as the weight 
of the piezo-actuator. 

Such a formulation results in a nonlinear optimization 
problem in 3 design variables t-j, v and i which is subjected to 
seven inequality constraints. The first inequality constraint 
guards against the application of a high voltage across the 
actuator that may result in its depolarization. The second 
inequality constraint limits the stresses in the piezo-film 
within the acceptable safe limits defined by the allowable stress 
a*. The last two inequalities limit the linear and angular static 
distortions to set permissible operational limits. 

A penalty function algorithm is used to transform the 
constrained minimization problem into an unconstrained 
optimization which is solved using a Powell's Conjugate Direction 
method [171. 

Figure (4) shows a block diagram of the optimization 


algorithm 



Input 

Beam geometrical and physical 
properties 


Input 

Initial location i, thickness t^ 
and voltage v of actuator 


Input 

Constraints on voltage, deflection, 
thickness and stresses 


Call the analysis algori 
compute Ur 


thm to 
) 


Check constraints and 
penalty function 

P= (Y i 2 +0 1 2 )+R*g 2 (x) 


Call Powell's Conjugate 
optimization algorithm 


/ are \ 

:onstraints satisfied 
and 

^stopping criterion^ 
\ met / 


Figure (4) - Flowchart of the optimization algorithm 






NUMERICAL EXAMPLES 


The developed optimization procedure is used to optimally 
select, place and control the action of a single piezo-electric 
actuator bonded to a three-element straight steel cantilever 
beam, shown in Figure (5)> that is 0.0125m wide, 0.0021m thick 
and 0.15m long. The effect of physical properties of the piezo- 
actuator and the bonding layer on the optimal configuration of 
the system is presented. 

The physical properties of the considered piezo-actuators 
and a typical bonding material are given in Tables (1) and (2) 
respectively . 

A. Performance With PZT Actuator 

m neglectin g the effect of bonding layer 

Figure (6) displays a 2-D plot of the design space in the 
t-|— v plane of the 3-D problem when one PZT piezo-actuator is 
placed at element 3 near the fixed end of the beam. This plot 
will give an understanding of the interaction between the 
constraints and the design objective. Plotted on this plane are 
the objective function contours and the constraints along with a 
search path to the optimum from a selected starting point. For 
such a case, the objective function assumes a minimum of 



Table (1) - Properties of piezo-electric actuators 


Actuator type 

Ceramic [Ref. 8] 

1 

Polymeric [Ref. 18] 

Material 

PZT 

1 

Kynar (PVF2) 

Charge coefficient d 
10” 12 (m/v) 

123 

1 

23 

Young's modulus 
(G N/m 2 ) 

139 

1 

1 

2 

Max. voltage v* 
(Mv/m) 

1 

1 

1 

30 

Max. tensile strength 
(M N/m 2 ) 

45 

1 

1 

33-55 

Density 

7500 

1 

1780 

(kg/nw) 


1 



Table (2) - Properties of a typical bonding layer [Ref. 10] 


Bonding Layer 

| Density (Kg/m^) | 

Young’s Modulus (MN/m^) 

Eastman 910^ 

1 1050 | 

1780.0 


1 Cyanoacrylate adhesive made by eastman Kodak Co., Rochester, 
NY 
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2.74x10“ 12 when a 0.08mm thick PZT actuator is placed at element 
3 and is energized by 80.00 volts. It can be seen that the 
optimum point lies, for the considered case and constraints, on 
the polarization constraint. 

Now to gain an insight into the effect that third dimension, 
i.e. the actuator location, has on the deflection characteristics 
of the beam let us consider the profiles of the different elastic 
lines obtained by varying the location of the actuator. The 
linear and angular elastic deformation are shown in Figures (7-a) 
and (7-b) respectively for both the uncontrolled and controlled 
beams. It is evident from the figures that placing an actuator at 
element 1 produced the least compensation for the external loads 
and its effects are limited to the first element only. However, 
placing the actuator at the third element, near the fixed end, 
results in considerable improvement that is not confined to 
element 3 is manifested clearly all over the three elements. 

(ii) Considering the ef fect of bonding layer. 

When the elastic and inertial properties of the bonding 
layer are accounted for, then the optimum values of the voltage, 
thickness of PZT actuator and the objective function are as shown 
in Figure (8) as function of the thickness of the bonding layer. 
The displayed results are obtained with Eastman 910 as the 
bonding layer. The figure indicates that increasing the thickness 
of the bonding layer results in reducing the optimum values of 





Applied Voltage (Volts) 



both the applied voltage (v) and the actuator thickness (t^) in a 
monotonic fashion. This can be attributed to the stiffening 
effect produced by the addition of bonding layer to the actuator- 
beam system as measured by the augmented flexible rigidity, 
computed from equation (8), and the element stiffness matrix 
given by equation (9). Such an effect results in reducing the 
deflections of the uncontrolled system. Accordingly, smaller 
voltages and actuator thickness would be needed to compensate for 
these reduced deflections. The figure indicates, however, that 
the effect of bonding layer on the optimal parameters v and t-j 
starts to level off beyond certain bonding layer thickness t 2 of 
0.00038m. This is due to the fact that the increased flexural 
rigidity of the composite is accompanied by a corresponding 
increase in the piezo-electric moment as given by equation 
(7). Therefore, these two factors balance each other to obtain a 
minimum objective function. 

B. Performance With Kynar Actuator 

Figure (9) shows a projection of the feasible design space 
in the v-t-j plane when one Kynar piezo-actuator is placed on 
element 3 without bond. 

The figure indicates that the objective function assumes a 
minimum of 1.24x10“9 when a 0.00244mm thick Kynar actuator is 
bonded to element 3 and is powered by 300 volts. 



OPTIMUM 
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Figure (9) - Iso-objectlve function contours with one Kynar 

actuator placed at element 3 neglecting effect of 
bond (V*s10 6 v/m, a* s45MN/m 2 , t^ # sO .0084m, 
v ma v=300v, W0.°0005m and ^sO.OOOSrad.) 



Therefore, under the same conditions and constraint limit, 
it is seen that the use of kynar actuator results in pushing the 
optimum towards the maximum voltage constraint instead of the 
polarization constraint as obtained with PZT actuators. More 
importantly, the optimum design with the Kynar actuator requires 
a piezo film which is about 3 times thicker than the optimal PZT 
film. Furthermore, in spite of the higher voltage and thicker 
actuator, the attained design objective is almost three orders of 
magnitudes higher than that obtained with the PZT actuator. 

Accordingly, from the v, t-j and U points of view, the use of 
the PZT actuator is favored over the Kynar actuator. 


Similar treads, to the PZT case, are also observed when the 
effect of the elastic and inertial properties of the bonding 
layer are considered as shown in Figure (10). 


CONCLUSIONS 


This paper has presented an optimum design procedure for 
controlling the shape of flexible beams using piezo-electric 
actuators. The procedure is based on accounting for the effect of 
the physical parameters of the actuator and its bonding layer on 
the elastic and inertial characteristics of the flexible beam. 


The presented procedure is utilized to select the voltage, 
thickness and location of piezo-actuators to minimize the 



Applied Voltage (10 Volts) 



structural distortion of flexible beams. 

The obtained results indicate that ceramic actuators are 
more effective than polymeric actuators in controlling the shape 
of a beam with lower voltage and thinner piezo films. Also, it is 
found that the addition of a bonding layer between the actuator 
and the beam results in stiffening the actuator-beam system and 
accordingly reduces the required applied voltage as well as 


actuator thickness 
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